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Growing evidence indicate supermassive black holes (SMBHs) in the mass range of Mbh^ 10 6 — 
1O 1O M lurking in central bulges of many galaxies [|lj ^j. Extensive observations reveal fairly 
tight power laws of Mbh versus the mean stellar velocity dispersion a of the host bulge [fl Q, 
The dynamic evolution of a bulge and the formation of a central SMBH should be physically 
linked by various observational clues. In this contribution, we reproduce the empirical Mbh — c 
power laws based on a self-similar general polytropic quasi-static bulge evolution [|(| 0] and a 
sensible criterion of forming a SMBH surrounding the central density singularity of a general 
singular polytropic sphere (SPS) [^||. Other properties of host bulges and central SMBHs are also 
examined. Based on our model, we discuss the intrinsic scatter of the Mbh — G relation and a 
scenario for the evolution of SMBHs in different host bulges. 
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Supermassive Black Holes and Galactic Host Bulges 

SMBHs form at the centres of elliptical and spiral galaxies JT0|, [T], p|]. Observationally, SMBH 
masses Mbh correlate with various properties of spiral galaxy bulges or elliptical galaxies, includ- 
ing bulge luminosities [jlj |l2|, 11], bulge masses Mb u i ge [12, 11, 13], galaxy light concentrations 

, inner core radii [|l(]], bulge gravitational 
I, Hi. These correlations 



[]14j], the Sersic index of surface brightness profile 
binding energies [|l7|] and mean stellar velocity dispersions a [18,19, 



strongly suggest a dynamical link between SMBHs and their host bulges [21 



Among these relations, Mbh and a correlate tightly in power laws with an intrinsic scatter 



< 0.3 dex [35]. This relation was studied [24, 23, 22fl before systematic observations [18, 19] and 



the emphasis was on outflow effects of galaxies. The idea was further elaborated [|25|]. A model 
of singular isothermal sphere (SIS) with rotation [26] was proposed for the Mbh — o relation. 



This relation was also explored in a semi-analytic model [27] with starbursts while SMBHs being 
formed and fueled during major mergers. Accretions of collisional dark matter onto SMBHs may 
also give the Mbh — <7 relation [28] (see also ref. [29] for a review). There are also numerical 



simulations to model feedbacks from SMBHs and stars on host galaxies. 

There are two empirical types of bulges, namely classical bulges (spiral galaxies with clas- 
sical bulges or elliptical galaxies) and pseudobulges [pC|]. While SMBHs in classical bulges are 
formed after major mergers, pseudobulges do not show apparent merger signatures. Interestingly, 
pseudobulges also manifest a Mbh — power law yet with a different exponent [Q]. 

The self-similar quasi-static solutions take the singular polytropic spheres (the static polytropic 
solutions) as the leading terms. This kind of dynamic solutions has been applied to study the 
compact stars in a single fluid model galaxy clusters in a two-fluid model [0] and the so-called 
"champagne flows" in H II regions In the following, we shall take the quasi-static solutions in 
a single fluid model to describe host bulges with central SMBHs after long time evolution [§]. 



A Self- Similar Dynamic Model for Mbh — o Power Laws 

For the dynamic evolution of a galactic bulge, we adopt a few assumptions. First, we treat 
the stellar bulge as a spherical polyuopic fluid as the typical age ~ 10 9 yr of galactic bulges is 



long [[31], |32j] that they are continuously adjusted. Stellar velocity dispersions produce an effective 



pressure P against the self-gravity as in the Jeans equation [|33p- Secondly, the total mass of the 
interstellar medium in a galaxy [g2j is ~ 10 7 - 10 8 M e only 1(T 2 ~ 10" 3 of the total bulge mass. 



Although gas densities in broad and narrow line regions of AGNs are high, the filling factor [ |34| ] is 
small (~ 10~ 3 ) and the gas there may be regarded as condensed clouds. Thus gas is merged into 
our stellar fluid. Thirdly, the diameter of broad line regions of AGNs [34] is only ~ 0.1 pc and 



the disk around a SMBH is even smaller while a galactic bulge size is > 1 kpc. We thus ignore 
small-scale structures around the central SMBH of a spherical bulge. Finally, as the rotation curves 
of galaxies show [[33|], the effect of dark matter halo in the innermost region (around several kpcs ) 
of a galaxy can be neglected. So we do not include the dark matter in our model when we discuss 
the dynamic characters of the bulges. 

Hydrodynamic equations of our model are conservations of mass, radial momentum [Q] and 
'specific entropy' along streamlines [^, ^8] [J. As bulk flow of stellar fluid is slow, we invoke the 
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novel self-similar quasi-static solutions [7|] to model the bulge evolution. We use a self-similar 
transformation [||, |8|, 0] to solve general polytropic fluid equations, namely 

r = K l ' 2 t n x, p = a(x)/(4xGt 2 ) , u = K l ' 2 t n - l v{x) , 
P = Kt 2n - 4 p(x)/{4nG) , M = K 3 / 2 t 3 "- 2 m(x)/[(3n-2)G], (1) 

where r is the radius and t is time; x is the independent dimensionless similarity variable while K 
and n are two scaling indices; G = 6.67 x 10~ 8 g _1 cm 3 s -2 is the gravity constant; p(r, ?) is the 
mass density and a(x) is the reduced mass density; u(r, t) is the radial flow speed and v(x) is the 
reduced flow speed; P(r, t) is the effective pressure and j8(x) is the reduced pressure; M(r, t) is 
the enclosed mass and m(x) is the reduced enclosed mass; reduced variables a, /3, v and m are 
functions of x only. We require n> 2/3 for a positive mass ||6[ [7|, 38]. 



By self-similar transformation ([]]), we readily construct self-similar quasi-static solutions from 
general polytropic fluid equations, taking the static singular polytropic sphere (SPS) as the leading 
term. Properties of such asymptotic solutions to the leading order ^ ^ are summarized below. 
Both initial and eventual mass density profiles scale as ~ r ^ 2 /"; accordingly, the bulge enclosed 
mass profile is M <=c r (3n-2)/«. f or e jt ne r x q+ or x — > +°°, the reduced velocity v — > 0, which 
means at a time f, for either r — > + or r — > +oo the flow speed m — > 0, or at a radius r, when t is 
short or long enough, the radial flow speed u — > 0. Our model [||] describes a self-similar bulge 
evolution towards a nearly static configuration after a long time lapse, appropriate for galactic 
bulges at present epoch. 

As the effective pressure P results from stellar velocity dispersions in the bulge, we readily 
derive the mean velocity dispersion a in a bulge. By specific entropy conservation along stream- 
lines, we relate P with p and M and derive the P profile from our quasi-static solutions J4C| , 0]. We 
take the local stellar velocity dispersion as Of,(r, t) = (yP/p) 1 / 2 where 7 is the polytropic index 
of our stellar fluid. To compare with observations, we derive the spatially averaged stellar velocity 
dispersion a in the bulge. The bulge boundary is taken as the radius r c where p drops to a value 
p c indistinguishable from the environment. For a class of bulges with same n, p c is regarded as a 
constant for a class of environments. One can show that within r c 

a = [3n 1+<? / 2 y 1/2 /(4"- \)]{4KG Pc ) {l - n V 2 A 3nc >l A K l l 2 = BK^ 1 , (2) 

where q = 2(n + 7- 2) /(3n - 2) and A = {n 2 ^ /[2(2 - n)(3n - 2)]}- 1 /(»-3"<?/2) . 

A SMBH forms at the centre of a galactic bulge that evolves in a self-similar quasi-static 
manner. Such a SMBH was formed by the collapse of collections of stars and gas towards the bulge 
centre and grows rapidly by matter accretions at an earlier phase [[R]]. As the growth timescale 
for SMBHs is only ~ 10 5 yr, our quasi-static solutions describe the relatively quiescent phase of 
galactic bulges after the formation of central SMBHs as a longer history of a bulge evolution. The 
stellar fluid made up of stars and condensed gas clouds has a slow bulk flow speed towards the 
central SMBH, sustaining a reservoir of mass accretion for the circumnuclear torus and/or disk. 
We now introduce the criterion of forming a SMBH. A SMBH mass Mbh and its Schwarzschild 
radius r s are related by Mbh = r s c 2 / (2G) where c is the speed of light. As we have pointed out, 
the bulge enclosed mass is M« r (3«-2)/n_ jf at a rac jj us r ^ m becomes M = r s c 2 /(2G), then a 
SMBH forms. Only those quasi-static bulges with n < 1 can form central SMBHs (see Figure pi); 
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Figure 1: The criterion of forming a central SMBH in a self-similar quasi-static bulge evolution. The enclosed mass 
power-law isM« r l) - 33 7 w j m n _ q 751 (sgiy curve). Meanwhile, we draw a straight dashed line M s = rc 2 /(2G) for the 
mass of SMBH versus the Schwarzschild radius r. Here at r s = 7.2 x 10 7 km, the straight line intersects the enclosed 
mass curve for a 2.45 x 10 7 M ra SMBH. A self-similar general polytropic quasi-static solution of n < 1 can form a SMBH 
at the bulge centre by invoking this criterion. 

we further derive Mbh x r s x 

K l/{2-2n) and the 

power law below 

(\ n/(2—2ri) / \ 3n— 2/(2— 2n) . ,, , . 
=b) (l) (3) 

where «Sf depends on c, G, n, 7, p c , and the exponent 1/(1 — «) > 3 as 2/3 < « < 1 is required. 

While the Mbh — power law is very tight with intrinsic scatter < 0.3 dex for SMBHs and host 
bulges [35], such scatter is large enough for different fitting parameters [ |T8| , |19[ |5j|. By equation 



(H), we have a natural interpretation for intrinsic scatter in the observed Mbh — power law. In 
our model, all bulges with the same n lie on a straight line with the exponent 1/(1 —n) as Figure [| 
shows. For a fixed n, different bulges are represented by different K values in transformation ([!]), 
leading to different Mbh and c However, for bulges with different n values, they lie on different 
lines. For elliptical galaxies or bulges in spiral galaxies, they appear to eventually take the self- 
similar evolution described above with a certain n value. But pseudobulges may take on different 
n values. Observationally, we cannot determine a priori the specific n value for a bulge but simply 
attempt to fit all bulges with a single exponent, which then contributes in part to intrinsic scatter in 
the observed Mbh — c power law. 

To show this, we fit three published Mbh — power laws in Figure |[ The first one is 
M B h = 1-2 X lO 8 M (a/2OO km s -1 ) 3 75 given in ref. [|^] with our parameters {n, 7, p c } being 
{0.733, 1.327, 0.47 M Q pc~ 3 }; and the five points correspond to K = {0.8, 1, 2, 3, 4} x 10 23 
cgs unit. The second one is log(MBH/M Q ) = 8.13 +4.021og(a/200 km s _1 ) given in ref. [j|] with 
our parameters {n, 7, p c } being {0.7512, 1.33, 0.0122 M pc }; and the seven points corre- 
spond to K = {1, 2, 4, 6, 8, 10, 20} x 10 22 cgs unit. The third one is log(M BH /M ) = 8.28 + 
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Figure 2: Power-law Mbh — o relations according to our general polytropic model for quasi-static self- 
similar evolution. Three index values n = 0.7537, 0.7512, 0.7333 are adopted for three different kinds (solid 



[gj], dashed [g], dotted lines [|19||) of Mbh — tr power laws given by equation (3). Given n, we calculate 
the SMBH mass Mbh an d the mean stellar velocity dispersion a for a certain K value in transformation. A 
bulge has different cr for different K. Each bulge-SMBH system is represented by a point in this display. 
All systems of same n give a straight line, while systems of different n correspond to different straight lines 
here. 



4.061og(a/200 km s l ) given in ref. [Q] with our parameters {n, y, p c } being {0.7537, 1.332, 
0.00364 M pc~ 3 }; and the six points correspond to K = {0.6, 0.9, 1.5, 2.5, 3.5, 4.5} x 10 23 cgs 
unit. Clearly, to fit all these points in Figure |2| with a single power law, we would get a different 
result with higher intrinsic scatter. In the three fitting examples, bulge inflow speeds of stellar fluid 
are slow (~ 0.1 — 1 km s _1 ), a feature of our self-similar quasi-static solutions. Near the SMBH 
boundary r s , the inflow rest mass-energy flux falls in the range of 10 40 - 10 45 erg s 1 in these ex- 



amples, sufficient to supply the observed X-ray luminosities [p6|]. There can be outgoing accretion 
shocks around a SMBH in these flows. As the age of galactic bulges [31] is so long (~ 10 9 yr) that 
such shocks should have gone outside bulges. 

Besides the Mbh — relation, observations reveal Mbh ^bulge w ^ tn ^buige being the bulge 
mass [|l3|] and Mbh x with E g being the absolute value of the bulge gravitational binding 



energy [17]. Using our criterion of forming a SMBH and the bulge radius r c , we derive a power 



law between Mbh and Mb u i ge as Mbh x ^buke 3 "^- ^ or n = 0.75, our result leads to relations in 



ref. [ ]26[ ] but for a nonisothermal general SPS. The bulge gravitational binding energy, without 
contributions from dark matter halo and a disk [jl7|], is E g « f Q c GMpAnrdr. For self-similar 
quasi-static solutions, we obtain Mbh x E^ 5 5 "\ 

As another class of bulges, pseudobulges are believed to have formed without merging in 
contrast to classical bulges. Nonetheless, pseudobulges follow a Mbh — power law but with a 
different exponent Q g]. In ref. [g], this M B h - o relation is found to be log(M B H/M ) = 7.5 + 
4.51og(a/200 km s _1 ). This conclusion is natural according to our model in that pseudobulges 
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may take a different self-similar quasi-static evolution for a different n value. Due to their different 
formation history, they show a different Mbh — G power law as observed. For {n, y, p c } being 
{0.7778, 1.34, 0.000426 M pc~ 3 }, our model can also fit this power law. 



Conclusions and Discussion 



The tight Mbh — G power laws and other relations among the SMBH mass Mbh and known 
properties of host bulges strongly suggest coeval growths of SMBHs and galactic bulges [^, |2^, 
27]. In our model, while forming a Schwarzschild SMBH at the bulge centre (e.g., by collapse 
of gas and stars or by merging), the spherical general polytropic bulge evolves in a self-similar 
quasi-static phase for a long time. We then reproduce empirical Mbh — G power laws. Different 
energetic processes appear to give rise to different scaling index n values. 

Besides classical bulges and pseudobulges, there are also 'core' elliptical galaxies, thought to 
have formed by 'dry' mergers. A steeper Mbh — G relation exists in these galaxies as compared 
to that for classical bulges [16]. This can be accommodated in our scenario that all hosts of 
SMBHs may finally evolve into self-similar quasi-static phase with different scaling parameters 
(e.g., different n for the slope and different p c for the normalization of the M B h — G relation). We 
thus provide a unified framework to model the the relatively quiescent evolution phase of SMBH 
host bulges and SMBH masses. As the observed Mbh — G relation for classical bulges is tight, the 
elliptical galaxies and spiral galaxies appear to take close n values for merging processes. 

In our model, n is a key scaling index to determine the exponent of the Mbh — G power law. 
The smaller the value of n is, the steeper the profile of the density is and the smaller the index of 
the Mbh — G relation is. If we think the SMBHs are formed by collapse of stars and gas and a 
less steeper density distribution may provide a more effective mechanism to form SMBHs, then we 
conclude that for a certain value of velocity dispersions, the smaller the mass of the initially formed 
SMBH is, the smaller the value of n is. 
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